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varying external driving force on interface growth. First, we derive a relation 



X 



We present a theoretical and numerical investigation of the effect of a time- 



between the roughening exponents which comes from a generalized Galilean in- 
variance, showing how the critical dimension of the model is tunable with the 
external field. We further conjecture results for the exponents in two dimensions, 
and find consistency with data obtained through simulations of two models we 
expect to be in the same universality class. Finally, we discuss how our results 
can be investigated experimentally. 
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Open systems pose among the most challenging problems in condensed-matter 
physics. The transient dynamics and rich steady-state behavior of systems rang- 
ing from dendrites to convective fluid cells provide a unique testing ground for 
the ideas and methods of nonequilibrium physics. Many such problems involve 
pattern formation where interfaces are present. A particularly well-defined prob- 
lem of this kind is the roughening of a growing interface A great deal of 
theoretical progress has been made on this problem by analytic and numerical 
methods, following the now-classic work of Kardar, Parisi, and Zhang (KPZ) 
f2j. They proposed a nonlinear differential equation to model the height h of a 
growing interface driven by an external flux of particles: 



where v and A are constants, and r\ is a random noise which is assumed to satisfy 
Gaussian statistics with (i](x,t)ri(x' ,t')) = 2D5 d ~ 1 (x — x')5{t — t'), where D is a 
constant, and the brackets denote a statistical average. The vector x determines 
positions in a (d — l)-dimensional plane of a full space r — (x,y), while t is 
time. The term proportional to A describes lateral growth of the interface, and 
plays an important role in roughening. This nonlinearity also implies the model 
has no free energy. 

Salient features of this system are the average velocity d(h)/dt ~ const., and 



the interface width W = \jh 2 — (h) 2 , where the bar denotes spatial average. The 

width obeys 



for late times and large length scales, where L is the linear size of the growing 
substrate, and / is a scaling function. In any dimension, Eq. (1) presents an 
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invariance, commonly referred to as "Galilean", providing an exact relationship 
between \ an d z. In addition, for d — 2, a fluctuation-dissipation theorem allows 
one to calculate the interface exponents, x — 1/2 and z = 3/2 Q. These are con- 
sistent with numerical simulations and differ from the trivial exponents obtained 
for A = 0: x = (3 — d)/2 and z = 2 ||. A renormalization-group analysis shows 
that the critical dimension of the model is d c = 3, but that it involves an infra-red 
unstable fixed point. Thus d c is thought to be a lower critical dimension for a 
transition which occurs in d > 3 between a strong coupling and a weak coupling 
behavior. 

Many justifications of this equation have been made in different contexts, par- 
ticularly for sputtered growth of epitaxial layers. Here we focus on the interface 
of a stable phase growing at the expense of a metastable phase . One example 
of such a system is the full field model of critical dynamics, called model A 0. 
There a nonconserved order parameter field ip(r,t) is the only dynamical mode. 
In equilibrium, at low enough temperatures, a stable interface exists between co- 
existing uniform phases ip « ±1. If a system is prepared in this fashion and a 
constant external field H which is conjugate to ip is applied, one of the two phases 
becomes metastable, and the interface advances through it. It can be shown that 
the dynamics of the interface separating those phases is described by the KPZ 
equation, with A oc H || . Systems such as the surfaces of growing dendrites also 
correspond to the KPZ equation if one takes the appropriate limit of model C, 
where in addition to ip, there is a conserved field coupled asymmetrically to the 
order parameter. In this case, the diffusion length of the conserved field must be 
much larger than the system size. 

The purpose of this Letter is to study the interesting possibility of a time 
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varying external field H(t) in model A. The mapping giving the KPZ equation 
still follows, but now with 

X(t) oc H{t). (3) 

Since H is an external field its time dependence can be prescribed in any fashion. 
Most importantly, we expect it to be experimentally tunable, as we shall discuss 
below. 

Assuming now a time- dependent X(t) in the KPZ equation ([]]), there are a 
number of interesting properties, three of which we now describe. Our main 
working hypothesis will be that the scaling relation (0) still remains valid. 

1. The equation has an invariance under a generalized "Galilean" transforma- 
tion 

x -> x - e I A(t') dt', and h -> h + e- x, (4) 

J const. 

where e is a constant vector of infinitesimal amplitude. This implies that 
X(t) is not renormalized under the application of the renormalization group, 
since it is directly involved in an exact invariance of the equation. Here we 
consider a power-law dependence of X(t) on time as: 

Kt) = p (5) 

where Ao is a constant, and report results for a > (we have also studied 
a < 0). The invariance above implies the novel relations: 

X + z = 2 + za, and d c = 3 — 4a, (6) 

where the tunable critical dimension is obtained from the trivial A = 
exponents. If the scaling hypothesis is valid, these are exact results. We 
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note that in the asymptotic regime, d{h(t))/dt ~ t a . This model has been 
studied previously by Lipowsky || in a different context. 

2. The importance of the result (6) lies in that a nonzero a could be sufficient 
to move to lower dimensions the phase transition which the usual KPZ 
equation exhibits only for d > 3. This will make it experimentally accessi- 
ble. For such a transition, at small Ao we should have the trivial exponents, 
which satisfy (2 — x)/z = min((d+l)/4, 1) , while beyond some larger value 
of Ao, (§) should hold, implying (2 — x)j z = 1 — a. In principle, it should 
also be possible to investigate this transition analytically, but we have not 
been able to do so. We intend to numerically study the possibility of this 
transition in the future. 

3. The nontrivial requirement of Galilean invariance implies that exponents 
are tunable in the presence of a time-varying field. In d = 2, a fluctuation- 
dissipation relation is satisfied both for A = and A = const, (i.e., for 
a = +oo and a = 0) giving \ — 1/2- We therefore conjecture that the 
value of x remains 1/2 for any intermediate A oc l/t a , with a > 0. This 
hypothesis and Eq.(^j) fix the values of the scaling exponents: 

x = I, and , = -» (7) 
A 2' 2(1 -a) V ; 

for d = 2, with < a < 1/4. For a > 1/4 one is above the critical dimension 

so that, at least for small Ao, we expect the ideal interface results with z — 2. 

It then follows that one could have a phase transition in two dimensions 

for a > 1/4, if d c remains a lower critical dimension. Of course, the fixed 

point could become infra-red stable for some a in < a < 1/4, so that 

there would be no transition in d = 2. In d = 3, the situation may be 



even more complicated. If x were to take the ideal interface value of for 
a > 0, with z determined from Eq. (|6]), the scaling exponents would jump 
discontinuously for a > 0. This question will be studied in more detail in 
future work. 

We have checked our hypothesis numerically with two models for d = 2. First, 
we have directly integrated the KPZ equation with a time-dependent A. Second, 



we have simulated the restricted solid-on-solid growth (GRSOS) model ]T0|] using 
a time- varying growth rate. 

The details of our numerical work are as follows. The KPZ equation, dis- 
cretized as in || was integrated by the Euler method. Systems sizes ranging from 
50 to 50000 lattice sites were studied. The time step was generally taken as 0.02, 
and some smaller values were considered to check accuracy. Typical values for the 
rest of the parameters were v = 3.5, Ao = 50 and D = 0.01, and several values of 
a were considered. For the GRSOS model, systems of sizes up to L = 50000 were 
simulated also. The algebraic decay of A was achieved by incorporating time- 
dependent condensation and evaporation rates, whose difference is proportional 
to A. For both models, the known limits for a = and a = oo were checked. 

Direct determinations of x were obtained by running systems of several sizes 
L until saturation, and then using the scaling of W(L,t = oo) with L. For both 
models, we obtained x ~ 0.5 for all values of a studied, as shown in Fig. 1. Our 
conjectures are then most simply tested through the exponent j3 — %j 'z, where 
W ~ t 13 if t « L z . For d = 2 we expect 

^ (8) 
for < a < 1/4 and (3 = 1/4 for a > 1/4. To obtain /3, we performed extensive 



calculations of the time-dependence of the width. Our results are summarized in 
Fig. 2. Although it becomes increasingly difficult to extract f3 around a = 1/4, 
our results for various values of a are in complete agreement with Eq. (|8|). 

A somewhat different way of checking our prediction for the exponents can be 
obtained from the crossover scaling form of W(L = oo,t) proposed in H], which 
for d = 2 reads: 

W(L = oo,t)~t$g(t\+) , (9) 
with g(x) ~ x 1//12 for small x. Theoretical arguments and computer simulations 



of discrete growth models [11] have so far given = 4, while direct simulations 



of Eq. ([[]) seem to give |6|, [12] = 3. If we assume that (|9]) remains valid for 
X(t) ~ t~ a , we find W ~ t( 4 -^)/ 12 , if a < 1/0. For = 4 this agrees with our f3 
in d = 2. In fact, we have checked the scaling form @, with a time-dependent A 
for both the GRSOS model and the KPZ equation, and obtained = 4 for both 
models. This lends strong support to our direct determination of (3 (see Fig. 2). 
To obtain this result for the KPZ equation, it is crucial to realize that the lattice 
constant of the mesh in which we solve the equation must be smaller than the 
length scale i> 3 /\ 2 D. For coarser discrete meshes, as used in the previous works 
|| |i~2| , we observe a crossover to = 3. Details of these results will be published 



elsewhere. 

Experimentally, the phenomena described here can be probed with metastable 
systems such as a solid growing into a supercooled melt or into a supersaturated 
solution. The degree of metastability, which is proportional to A, can be varied by, 
for example, progressively decreasing the degree of undercooling. Similarly, in a 
sputtered growth experiment, our results describe a system as the sputtering rate 
is tuned to fall off algebraically in time. Thus we expect it to be experimentally 
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feasible to probe the three properties we have discussed: The generalized relation 
of Galilean invariance, the possibility of a phase transition at physical dimensions, 
and the tunable roughening exponents. 

This work was supported by the Natural Sciences and Engineering Research 
Council of Canada, and le Fonds pour la Formation de Chercheurs et I' Aide a la 
Recherche de la Province de Quebec. E.H-G acknowledges support from Comision 
Interministerial de Ciencia y Tecnologia (Spain). T.A-N. was also supported by 
the Academy of Finland, and acknowledges the hospitality of McGill University 
during his visit there. 



8 



References 

[1] J. Krug and H. Spohn, in Solids Far From Equilibrium: Growth, Morphology 
and Defects, edited by C. Godriche (Cambridge Univ. Press, New York, 
1990). 

[2] M. Kardar, G. Parisi and Y. C. Zhang, Phys. Rev. Lett. 56, 889 (1986). 

[3] Here, d is the dimension of the full space. A common notation is to write 
d — (d — 1) + 1. In our notation 1 + 1=2, 2 + 1 = 3, and so on. 

[4] D. Forster, D.R. Nelson, and M.J. Stephen, Phys. Rev. A 16, 732 (1977). 

[5] S.F. Edwards and D.R. Wilkinson, Proc. R. Soc. London, Ser. A 381, 17 
(1982). 

[6] B. Grossmann, H. Guo, and M. Grant, Phys. Rev. A 43, 1727 (1991). 

[7] PC. Hohenberg and B.I. Halperin, Rev. Mod. Phys. 49, 435 (1977). 

[8] The mapping follows on large time scales which are, nevertheless, sufficiently 
small that nucleation does not take place in the metastable state. 

[9] R. Lipowsky, in Random Fluctuations and Pattern Growth: Experiments 
and Models, edited by H.E. Stanley and N. Ostrowsky (Kluwer, Dordrecht, 
1988). This model was studied in this reference to investigate wetting of a 
substrate. That involves particular values for the exponent a, for which A 
turns out to be an irrelevant variable. 

[10] J. M. Kim and J. M. Kosterlitz, Phys. Rev. Lett. 62, 2289 (1989); J. M. 
Kim, J. M. Kosterlitz, and T. Ala-Nissila, J. Phys. A 24, 5569 (1991). 

9 



[11] J. Krug, P. Meakin and T. Halpin-Healy, Phys. Rev. A 45, 638 (1992); J.G. 
Amar and F. Family, ibid., 5378 (1992). 

[12] B.M. Forrest and R. Toral, J. Stat. Phys. , in press (1992) 



10 



Figure Captions 

1 . The exponent x as a function of a from simulations of the KPZ equation (di- 
amonds) and the GRSOS model (asterisks). KPZ data come from averages 
over 100 realizations of systems of sizes L = 50, 100, 200, and 400. GRSOS 
data have been averaged over 10 7 — 10 8 Monte Carlo steps for L = 100, 200, 
500 and 1000. The error bars are purely statistical; no systematic finite size 
errors have been included. 

2. The exponent (3 as a function of a. Meaning of the symbols as in Fig. 
1. Solid line is the prediction of Eq. (§). Data come from simulations of 
systems of size L = 50 000 averaged over 30 (KPZ) and 30-60 (GRSOS) 
independent runs. 
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